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ABSTRACT 

We present the energy-dependent power spectral density (PSD) and cross-spectral properties of Mkn 
766, obtained from combining data obtained during an XMM-Newton observation spanning six revolu- 
tions in 2005 with data obtained from an XMM-Newton long-look in 2001. The PSD shapes and rms-flux 
relations are found to be consistent between the 2001 and 2005 observations, suggesting the 2005 ob- 
servation is simply a low-flux extension of the 2001 observation and permitting us to combine the two 
data sets. The resulting PSD has the highest temporal frequency resolution for any AGN PSD measured 
to date. Applying a broken power-law model yields break frequencies which increase in temporal fre- 
quency with photon energy. Obtaining a good fit when assuming energy-independent break frequencies 
requires the presence of a Lorentzian at 4.6 ± 0.4 x 10~ 4 Hz whose strength increases with photon en- 
ergy, a behavior seen in black hole X-ray binaries. The cross-spectral properties are measured; temporal 
frequency-dependent soft-to-hard time lags are detected in this object for the first time. Cross-spectral 
results are consistent with those for other accreting black hole systems. The results are discussed in the 
context of several variability models, including those based on inwardly-propagating viscosit}^ variations 
in the accretion disk. 

Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies — galaxies: individual (Mkn 
766) 


1. INTRODUCTION 

Seyfert active galactic nuclei (AGNs) and stellar-mass 

black hole X-ray binary systems (XRBs) both exhibit 

rapid, aperiodic X-ray variability that likely originates 

in the innermost regions of these compact accreting ob- 
jects. The dominant X-ray radiation is generally thought 
to be inverse Comptonization of soft seed photons by a 
hot corona (e.g., Shapiro et al. 1976, Sunyaev & Titarchuk 

1980), though the exact geometry is uncertain and numer- 
ous configurations have been invoked (e.g., Zdziarski et 
al. 2003). The X-ray variability can be characterized by 
fluctuation power spectral density functions (PSDs) which 
reveal the ’’red- noise” nature of the variability at rela- 
tively high temporal frequencies in both AGNs and XRBs. 
Modeling of the broadband PSDs of XRBs usually utilizes 
some combination of one, two, or more Lorentzians, which 
tend to dominate in the so-called low/hard energy spectral 
state, plus a broadband noise component characterized as 
P(f) oc / _1 , which tends to dominates in the high/soft 
state (e.g., in Cyg X-l; see Axelsson et al. 2005). Quasi- 
periodic oscillations, or QPOs, are detected in some XRBs, 
usually those in the so-called Very High spectral state (e.g., 


McClintock & Remillard 2006). 

Broadband PSDs have also been constructed for AGNs, 
allowing modeling of the overall PSD shape. Specifically, 
’breaks’ in the power-law PSD models on time scales of 
a few days or less have been identified for roughly a 
dozen AGN (e.g., Uttley, M c Hardy & Papadakis 2002; 
Markowitz et al. 2003; Vaughan, Fabian & Nandra 2003b, 
hereafter VFN03; M°Hardy et al. 2004, 2005), with power- 
law slopes < —2 above the break and ~ —1 below it. 
No evidence for QPOs has been found in AGN PSDs yet. 
One reason for this may that the insufficient sampling of 
the observed AGN light curves so far does not allow one 
to achieve a high enough temporal frequency resolution 
(Vaughan h Uttley 2005). Nonetheless, the overall simi- 
larity in Seyferts’ and XRBs’ broadband PSD shapes and 
scaling of PSD break time scale with black hole mass (e.g., 
Markowitz et al. 2003, M c Hardy et al. 2004), as well as 
the presence of a linear relation between flux and absolute 
rms variability (e.g., Uttley & M c Hardy 2001, Edelson et 
al. 2002), support the notion of similar X-ray variabilhy 
mechanisms being present in both classes of accreting com- 
pact objects. That is, Seyferts and XRBs may simply be 
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scaled versions of each other in black hole mass and X-ray 
variability time scale. 

Further support for this picture comes from observations 
of the PSD slope above the break to flatten with increasing 
energy. Previous studies of both XRBs (e.g., Nowak et al. 
1999a, Lin et al. 2000) and Seyferts (Nandra & Papadakis 
2001, VFN03, M c Hardy et al. 2004, and Markowitz 2005; 
see also Leighly 2004) have suggested that the PSD slope 
flattens as photon energy increases. There has been no ev- 
idence yet for break frequency to change significantly with 
photon energy in either AGN or XRBs during a single ob- 
servation. 

The observed cross-spectral properties, namely the co- 
herence function (Vaughan & Nowak 1997) and time lag 
spectrum, are also similar between Seyferts and XRBs. 
The coherence is generally seen to be close to unity 
over a wide range of temporal frequencies in both AGNs 
and XRBs, usually dropping significantly below unity 
at frequencies above the PSD break (e.g., Nowak et al. 
1999a, 1999b in XRBs; VFN03, M c Hardy et al. 2004 and 
Markowitz 2005 in AGN). Relatively harder X-rays are 
seen to generally lag those at softer X-ray energies, with 
the lag between bands increasing with increasing energy 
separation. However, lags are are also observed to increase 
with decreasing temporal frequency (e.g., Miyamoto & Ki- 
tamoto 1989, Miyamoto et al. 1991, Nowak et al. 1999a in 
XRBs), with time lags usually < 1 — 10% of the Fourier 
time scale. In the low-hard state of Cyg X-l, the time lags 
increase in a step-like fashion (e.g., Nowak 2000), which 
can be explained if the lag at each step is associated with 
individual Lorentzian components in the broadband PSD. 
Obtaining time lags via the cross-spectrum requires ex- 
tremely high quality data with sufficient sampling and 
number of points. Results for AGN so far have yielded 
time lag spectra with many fewer points than in XRBs, 
but also showing the general increase in time lags as tem- 
poral frequency decreases (Papadakis, Nandra & Kazanas 
2001, VFN03, M c Hardy et al. 2004, Markowitz 2005, and 
Arevalo et al. 2006a). 

Mkn 766 is a well-studied Narrow-Line Seyfert 1 
(NLSyl) , and thus belongs to a class of objects in which 
large-amplitude rapid X-ray variability is commonly ob- 
served. Timing properties of Mkn 766 have been studied 
previously by Vaughan <fc Fabian (2003; hereafter VF03), 
who measured the broadband PSD, finding a break from 
steep PSD power-law slope of < —2.5 to ~ —1 near 
5 x 10~ 4 Hz. VF03 also observed the slope above the 
break to flatten from ~ —2.7 at soft X-rays to ~ -2.2 at 
hard X-rays; however, they saw no obvious trends in break 
frequency with energy. VF03 were the first to measure the 
coherence in Mkn 766. In addition, Vaughan et al. (2003a) 
observed the linear rms-flux relation in Mkn 766. 

In this paper, we present the timing properties of Mkn 
766 using the enormous wealth of data obtained from un- 
interrupted XMM-Newton observations in 2005 and 2001. 
We find that the variability properties of the source during 
both observations were identical. We therefore combined 
the 2001 and 2005 light curves to construct power- and 
cross-spectra with the highest frequency resolution in log 
space seen an in AGN PSD so far. The rest of this paper 
is organized as follows: §2 describes the data reduction. 
§3 demonstrates similarities in the variability properties of 


the 2001 and 2005 observations. §4 describes measurement 
of, and model fits to, the energy-dependent PSD made 
from combining all the 2001 and 2005 data. §5 describes 
the cross-spectral properties. The results are discussed in 
§6 in the context of phenomenological variability models. 
Finally, a brief summary is given in §7. 

In one of two companion timing papers, Papadakis et al. 
(in prep.; hereafter Paper II) combine the 2001 and 2005 
XMM-Newton light curves with long-term light curves ob- 
tained from RXTE monitoring. That paper concentrates 
on study of the low-frequency PSD shape of Mkn 766 to 
probe similarities to the broadband PSD shapes of XRBs. 
In the other paper, Arevalo et al. (in prep.; Paper III) con- 
tinue the study of the energy-dependent variability prop- 
erties via Fourier-resolved spectroscopy. 

2. OBSERVATIONS AND DATA REDUCTION 

Mkn 766 was observed by XMM-Newton on 2005 May 23 
- 2005 June 3, over revolutions 999-1004. This paper uses 
data taken with the European Photon Imaging Camera 
(EPIC), which consists of one pn CCD back-illuminated 
array sensitive to 0.15-15 keV photons (Struder et al. 
2001), and two MOS CCD front-illuminated arrays sen- 
sitive to 0.15-12 keV photons (MOS1 and MOS2, Turner 
et al. 2001). Data from the pn were taken in Small Window 
mode; data from each MOS were taken in Large Window 
mode. The medium filter was used for all detectors. 

Light curves were extracted using SAS version 6.5 using 
standard extraction procedures. Data were selected using 
event patterns 0-4 for the pn and 0-12 for each MOS. Hot, 
flickering or bad pixels were excluded. Source light curves 
were extracted from circular regions of radius 40"; back- 
ground light curves were extracted from circles of identical 
size, centered ~ 3' away from the source. We checked for 
pile-up; there was no significant pile-up in the pn and at 
most a few percent pile-up in each MOS. 

The 2-12 keV pn light curve was analyzed for back- 
ground flares; any data where the pn background rate 
exceeded 0.2 ct s _1 were excluded. Most data excluded 
occurred at the beginnings or ends of revolutions (e.g., the 
final 18 ksec of revolution 999 data, the final 17 ksec of 
revolution 1003, and the final 13 ksec of revolution 1004). 
This screening also introduced a 10 ksec gap into the rev- 
olution 1003 light curve. Because such large gaps tend to 
reduce the statistical significance of variability parameters 
derived over the full duration, the light curve was split into 
two smaller light curves, hereafter designated 1003 A and 
1003B, for a total of seven virtually uninterrupted light 
curves. All other gaps introduced by background screen- 
ing were small (< a few ksec) and in the analyses below, 
fluxes during gaps were linearly interpolated from adjacent 
points. 

Data were extracted over the 0.2-12 keV (total), 0.2-0. 5 
keV (soft), 0. 5-1.1 keV (medium), and 1.1-12 keV (hard) 
bands (also referred to as T, S, M and H bands, respec- 
tively, below). The S, M and H band limits were cho- 
sen such that all three bands had approximately the same 
count rates. Consequently, the variability due to Poisson 
noise (see §4 and 5) was roughly the same level, mitigating 
any effect that the difference in count rates or Poisson noise 
between bands could have on analysis. Light curves from 
the three EPIC cameras were summed, using only data 
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taken when all three cameras were in operation simulta- 
neously; this yielded a total exposure time (summed over 
all light curves) of, 442 ksec. Given the wide bands used 
here, the differences between pn and MOS response shapes 
have a negligible effect; the average pn and MOS photon 
energies for each band differed only by small amounts. In 
addition, Mkn 766 is a soft-spectrum source, and the in- 
clusion of the MOS data compared to the pn data alone 
greatly improved signal/noise in the M and H bands. Light 
curves were binned to 60 s; they are shown in Figure 1. 
The summed T, S, M and H band light curves have mean 
count rates (after background subtraction) of 13.44, 4.39, 
4.39, and 4.93 ct s -1 , respectively. We use these summed 
light curves in all analyses below except the rms-flux rela- 
tion and F var spectra, where only pn data are used. 

Mkn 766 was also observed during revolution 265 from 
2001 May 20 - 2001 May 21. The pn and MOS2 were both 
in Small Window mode; the MOS1 was in Timing mode 
and is not considered further. Data were extracted from 
the pn and MOS2 using the same procedure and bands 
as for the 2005 data. No significant pile-up was found in 
the pn or MOS2. This reduction yielded a 107.1 ksec light 
curve summed over the pn and MOS2 The summed T, S, 
M and H band light curves have mean count rates (after- 
background subtraction) of 26.85, 9.18, 9.34, and 8.12 ct 
s _1 , respectively. The pn light curve is presented in VF03. 

3. BASIC VARIABILITY PROPERTIES 

To quantify variability amplitude as a function of en- 
ergy, we constructed F v ar spectra for the 2005 observation 
using the pn data and bin sizes of 6 ksec, and using the 
formulation of Vaughan et al. (2003a). Figure 2 shows the 
F V ar spectrum over the entire observation, and during the 
lowest- and highest-flux revolutions; the plots confirm rel- 
atively high levels of variability, as expected for a NLSyl. 
The strongest variability amplitudes are found near 1-2 
keV in both the 2001 (VF03) and 2005 F var spectra, a re- 
sult similar to e.g., VFN03 for MCG-6-30-15. The results 
corroborate the principal component analysis of Miller et 
al. (2006), who demonstrate that the variability of the con- 
tinuum component in the soft and hard bands is diluted 
by the presence of less-variable components, namely the 
Compton reflection spectrum in the hard band and by the 
soft excess and narrow emission lines in the soft band. 

Compact accreting black holes display “weakly non- 
stationary” behavior in that the mean and variance both 
show scatter over time, although both the underlying PSD 
shape and the expectation value of the fractional variabil- 
ity amplitude <F var > are expected to remain constant in 
time, e.g, over time scales shorter than decades to cen- 
turies in AGN. A constant <F var > in flux is equivalent to 
the linear rms-flux relation; Vaughan et al. (2003a) pre- 
viously established the presence of both a linear rms-flux 
relation and constant <F var > with flux in Mkn 766 using 
the 2001 data. 

We use the enormous wealth of data from the 2005 ob- 
servation to revisit this issue. The average 0.2-12 keV 
flux was 3 x 10 _n erg cm -2 s” 1 (e.g., Turner et al. 2006), 
about half that of the 2001 observation, 6 x 10 -11 erg cm -2 
s _1 (e.g., Pounds et al. 2003), allowing us to test if the 
rms-flux relation has persisted over a four year time scale 
and can extend to lower fluxes. We measured variabil- 


ity amplitudes in the T-band light curve over short time 
scales, using only pn light curves to facilitate compari- 
son between 2001 and 2005. We calculated the flux, F var , 
and the rms amplitude a rms every 2000 seconds (subtract- 
ing off the variance due to Poisson noise). Light curves 
with <10 points were discarded. Flux, rms amplitude, and 
F var points for all seven light curves from 2005 were com- 
bined; points were sorted according to flux, and binned in 
groups of 10. Results are shown in Figure 3. Both obser- 
vations are consistent with a linear <cr rms >-flux relation, 
as shown in top panel, and with <F var > being constant in 
flux, as shown in the bottom panel. Moreover, the 2001 
and 2005 relations are roughly consistent with each other, 
indicating that characteristics of the underlying variabil- 
ity process have not significantly changed since 2001, and 
the 2005 observation is merely a low-flux extension of the 
2001 observation. The positive x-intercepts of the best-fit 
<<7rms>-flux relations suggest the presence of a constant 
component whose magnitude has not changed significantly 
between the two observations; this qualitatively corrob- 
orates the principal component analysis of Miller et al. 
(2006). 

We performed another test for strong non-stationarity 
using the PSD shapes; as described in the Appendix, the 
PSD shapes from 2001 and 2005 are consistent with each 
other. Having established that the variability processes 
are consistent between the 2001 and 2005 observations, 
we combine the 2001 and 2005 light curves in all analy- 
ses below to maximize the temporal frequency resolution. 
The difference in mean count rates (e.g., summed over pn 
4- MOS2 in 2001 compared to pn + MOS1 -I- MOS2 in 
2005) is irrelevant since light curves are mean-subtracted 
and mean-normalized in all analyses below. 

Virtually identical results are obtained when using the 
2005 data only, though below we present results from 
the combined 2001 and 2005 binned PSD, which features 
~ 25% higher temporal frequency resolution in log space. 

4. THE ENERGY DEPENDENCE OF THE HIGH-FREQUENCY 

PSD 

The high-quality (uninterrupted and highly variable) 
light curves allows us to create a very high quality PSD 
for Mkn 766; in particular, the temporal frequency resolu- 
tion allows very tight constraints on break frequency. §4.1 
describes how the PSDs were constructed; §4.2 describes 
model fits to the PSD shapes. 

4.1. PSD Construction 

The construction of the PSDs is summarized here; for 
further details on PSD construction, see e.g., Uttley et 
al. 2002, Markowitz et al. 2003, VFN03, and references 
therein. Unlike Paper II, we do not include RXTE Pro- 
portional Counter Array (PCA) monitoring data, since the 
PCA bandpass does not go below 2 keV and this paper’s 
aim is to explore the PSD behavior across as wide an en- 
ergy range as possible. First, for each of the seven light 
curves from 2005 and for the 2001 light curve, the mean 
was subtracted. Periodograms were then calculated sep- 
arately for each of the eight light curves using a Discrete 
Fourier Transform (DFT; e.g., Oppenheim & Shafer 1975). 
The PSD normalization of van der Klis (1997) was used. 
All periodogram points from the eight light curves were 
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combined and sorted in temporal frequency. Following Pa- 
padakis and Lawrence (1993), the periodogram was binned 
in log space in groups of 25. This process yielded a PSD 
for each of the T, S, M and H bands, each with a usable 
temporal frequency range of 2.57 x 10 -5 — 8.27 x 10 -3 Hz. 
Power at temporal frequencies above ~3x lO -3 Hz were 
dominated by the white noise power due to Poisson noise; 
this observed level of power agreed well with predicted val- 
ues calculated as Pp sn = 2 (/i + where /jl and B are 

the exposure-weighted net source and background rates, 
respectively. The T-band PSD P(f) is plotted in Figure 
4a (in / — P(f) space). The S, M, and H-band PSDs are 
plotted in Figures 5a and 5b (both in / — fxP(f) space). 

4.2. PSD Model Fits 

We now discuss model fits to the broadband PSD shape. 
Monte Carlo simulations (e.g., Uttley et al. 2002) are fre- 
quently employed in PSD measurement when it is neces- 
sary to account for PSD measurement distortion effects 
and to assign errors to poorly-sampled PSD points. How- 
ever, for the present PSDs, we expect no aliasing since each 
of the individual light curves were continuously sampled 
and large gaps were excluded. The effects of red noise leak 
from temporal frequencies lower than those sampled here 
are likely small given the measured PSD shape, and at any 
rate are assumed not to vary significantly with photon en- 
ergy. Finally, because each PSD bin contains a sufficient 
number of unbinned periodogram points (25), PSD errors 
are normal, removing the need for Monte Carlo simula- 
tions to determine PSD errors. All model fitting was done 
in log space. We used least-squares fitting (e.g., Beving- 
ton 1969) to determine the PSD amplitude. Errors below 
correspond to A% 2 = 2.71. 

4.2.1. Unbroken and Broken Power-law Model Fits 

We first fit an unbroken power-law of the form P(/) = 
A(f/f 0 )-P + Pp sn where the normalization A is the PSD 
amplitude at the (arbitrary) frequency /o- /? is the power 
law slope; e.g., positive values of (3 denote red-noise PSD 
slopes. Ppsn is the constant level of power due to Poisson 
noise, kept fixed in all fits at the predicted values, which 
are calculated to be 0.16, 0.44, 0.48 and 0.44 Hz -1 for the 
T, S, M and H bands, respectively. Pp sn was added to the 
model as opposed to being subtracted from the data, to 
avoid the possibility of obtaining unphysical negative pow- 
ers. Values of (3 were tested in increments of 0.02. Results 
for the best-fit model are listed in Table 1; residuals to the 
fits are shown in Figures 4b (T-band) and 5c (S, M, and 
H bands). As per VF03, the fits are quite poor and signal 
the need for a break in the PSD model. 

We next tried fitting broken power-law models. In AGN 
PSDs, it has not always been clear if the breaks are sharp 
or instead follow a slow bend. We fit a sharply-broken 
power-law model of the form 

pm _ / A(f/f h )-T + Pp m , f<h 
\ A(f/f h )-e + Pp sn , f>f h 

Here, 7 and f3 are the power-law slopes below and above 
the break frequency / b ; A is the PSD amplitude at / b . We 
also tested a slowly-bending power-law model, 

P(/) = A/^(l + (/// b )^) + P Psn 


Here, 7 and /3 are the power-law slopes far below or above, 
respectively, the rollover frequency / b . A is the PSD nor- 
malization. Values of log(/ b ) were tested in increments of 
0.01; values of / 3 and 7 were tested in increments of 0.02. 
Results for the best-fit models with /?, 7, / b and A free are 
listed in Table 2; and the residuals are shown in Figures 
3c and 3d (T-band) and Figures 5d and 5e (S, M, and H 
bands). 

These fits are significantly better (all > 4.2<r confi- 
dence in an F-test) compared to the unbroken power-law 
model fits. Both sharply-broken and slowly-bending model 
shapes yield similar fits. There is no obvious trend in either 
/3 or 7 with photon energy in either case. We note that 7 is 
similar to the PSD power-law slope < 10 -3 ' 5 Hz obtained 
from fitting the XMM-Newton and long-term RXTE data, 
as shown in Paper II. For either shape, though, / b for the 
S band is significantly lower by ~ 0.2 — 0.5 in the log (de- 
pending on the model used) compared to / b for the M or 
H bands. 

We explored if it was significant to keep /3 and 7 thawed 
in the fits. We repeated the fits, with (3 and 7 frozen at 
the average value of the S, M and H-band fits. Results are 
listed in Table 2 and the best-fit models are plotted in Fig- 
ures 5a and 5b. The sum of S, M and H x 2 /dof (degrees 
of freedom) increased from 574.9/534 to 584.9/540 for the 
sharply-broken model, and from 577.6/534 to 588.3/540 
for the slowly-bending model. It was thus not significant 
at greater than the 87% confidence level (in an F-test) to 
thaw /? and 7. With (3 and 7 frozen, / b , which has typical 
errors of 0.05 in the log, is seen to increase with increas- 
ing photon energy (e.g., errors do not overlap at the 90% 
confidence level between the M and H bands). 

Refitting the PSDs with [3 , 7 and / b all kept frozen at 
their average values between the S, M and H bands yields 
a total (S+M+H) x 2 /dof of 604.2/543 (sharply-broken) 
and 615.9/543 (slowly-bending). Compared to the model 
fits with just f3 and 7 frozen, we find that it is significant 
at >99.99% confidence in an F-test (for either model) to 
keep / b thawed in the fits. 

We conclude that the data are not good enough to allow 
us distinguish, formally, between a sharp break or a slow 
rollover in the broken power-law models. However, the 
data are of high enough quality to conclusively demon- 
strate a rather unexpected result, one not observed pre- 
viously in the PSDs of either AGN or XRBs: the break 
frequency increases with photon energy. This increase in 
fb with photon energy holds regardless of whether we use 
sharply-broken or slowly-bending power law models and 
whether the power-law slopes in these models are frozen 
or thawed in the fits. Furthermore, the data are consis- 
tent with a broadband PSD shape that is identical for 
all bands in that the slopes above and below the break 
are energy-independent; this also has not been previously 
demonstrated. 

Depending on whether the sharply-broken or slowly- 
bending model is used, and assuming (3 and 7 are energy- 
independent, we find that from the S to M bands, / b in- 
creases by 0.08±0.08 or 0.05±0.09 in the log; from the M 
to H bands, / b increases by 0.13T0.07 or 0.16±0.06 in the 
log; and from the S to H bands, / b increases by 0.21±0.08 
or 0.21±0.08 in the log. In other words, the change in 
fb from the S to H bands is a 2.6cx effect. However, this 
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is a conservative estimate, since the errors in Table 2 are 
90%, not la, so the significance of this result is proba- 
bly larger than 2.6u. The average photon energy for each 
band, weighted for instrument count rate and exposure 
time, is 0.36, 0.75, and 3.9 keV for the S, M and H bands, 
respectively. The increase in fb between the S and H bands 
is a factor of ~ 1/5 times the ratio of the average photon 
energies. 

4.2.2. The Case for a Lorentzian peak in the PSD of 

Mkn 766 

Previous PSD papers have reported that /3 flattens as 
photon energy increases while fb is energy-independent. 
We fit the S, M and H bands PSDs with the same sharply- 
or slowly-bending power-law models, but with 7 and fb 
assumed to be identical for all bands. The best-fit results, 
calculated by minimizing the sum of x 2 from the S, M and 
H bands, are listed in Table 2. The sum of x 2 / do f for the 
S, M and H bands was 597.5/540 for the sharply-broken 
model and 599.1/540 for the slowly-bending model. The 
summed values of x 2 /d°f when fb was allowed to vary in 
photon energy (7 and (3 frozen) was 584.9/540 (sharply- 
broken) or 588.3/540 (slowly-bending). We conclude that 
a simple broken power-law model where fb does not vary 
with energy cannot accurately describe the S, M and El- 
band PSDs of Mkn 766. 

However, there is an alternate model that fulfills the. re- 
quirements of having fb energy-independent and providing 
a good fit to the S, M and H-band PSDs. This is a sharply- 
broken power-law model with the addition of a Lorentzian 
component whose strength increases in photon energy: 

pm _J Mf/hr^ + PL + Ppsn, f<fb 

~ \ A(f/f h y0 + P h + Ppsn: f>fb 

where 

R 2 hQ/n 

L /! + Q 2 (/-/l) 2 

fb is the Lorentzian centroid frequency. Q is the coher- 
ence or quality factor. R is the normalization factor of the 
Lorentzian; for high Q values, R is the fractional rms vari- 
ability amplitude of the Lorentzian. For values of Q < 2, 
the Lorentzian is very broad and cannot be easily distin- 
guished from band-limited or broadband noise, e.g., the 
broken power-law component. Lorentzians with higher Q 
values are called QPOs (e.g., van der Klis 2006). 

The addition of the Lorentzian component to the 
sharply-broken power-law model complicates the model- 
fitting process. For that reason, in order for the model- 
fitting results to be reasonable and of some use, we de- 
cided to make the following plausible assumptions regard- 
ing the parameter values of the Lorentzian component, 
based mainly on previous results from similar model fits to 
the XRBs PSDs: First, we are interested in fitting a peak- 
like component, and so we only consider a Lorentzian with 
Q > 3 in order to avoid fitting a component that resem- 
bles additional broadband noise. We cannot place tight 
constraints on Q but have found that values of ~ 4 — 5 
were the most plausible in terms of providing a good fit to 
the PSDs. For that reason, all model fits were performed 
with Q kept fixed at 5. Second, we forced /l to be the 
same in all three PSDs (i.e., we assumed that fb is intrin- 
sically energy-independent) and we further constrained its 


value to lie in the range 10~ 3 - 22 — 10~ 3 - 50 Hz; fb was tested 
in increments of 0.02 in the log. Finally, R was left as a 
free parameter, tested in increments of 0.5%. The break 
frequency fb of the broken power-law continuum compo- 
nent was not constrained to be equal for all bands, but 
was constrained to be < io -3 - 60 Hz. Finally, we left (3 as 
a free parameter but kept 7 fixed at —1.46 for simplicity. 

Table 3 lists the best-fit results obtained from finding the 
minimum value of x 2 summed over the S, M and H bands. 
The data, best fit models, and residuals are plotted in 
Figure 6. First, we note that both fb and the power-law 
slope above the break (3 are similar for all bands. Sec- 
ond, we find that the best-fit log(/L (Hz)) is — 3.34 ± 0.04 
(4.6 ± 0.4 x 10 -4 Hz). Finally, we note that R increases 
with photon energy: it is formally consistent with an up- 
per limit of only 1.5% in the S band, but is 5.5 ± 0.5% in 
the H band. This is also demonstrated by Figure 7, which 
shows derived confidence contours of R versus /l for each 
band separately. 

Using a slowly-bending model to fit the continuum 
yielded virtually identical values of x 2 > R , /l and iden- 
tical uncertainties on R and /l in each band; the best-fit 
value of fb was only 0.02 lower in the log compared to 
using the sharply-broken model. 

The temporal frequency resolution of this PSD, despite 
being the highest for any binned AGN PSD to date, is 
not high enough to break model degeneracies and simul- 
taneously constrain the curvature of the break and the 
profile of the Lorentzian component, especially since the 
Lorentzian is located so close to the break frequency. In 
addition, the broken power-law model and the model with 
the Lorentzian component both yield similarly good fits to 
the PSD, so we cannot exclude one or the other based on 
the quality of fit alone. 

Nonetheless, the PSD is of sufficient quality to detect 
energy-dependent behavior that has not been observed 
previously in an AGN. Either fb increases with photon en- 
ergy (in which case the broadband PSD shape is likely the 
same for all energies), or fb and /3 are energy-independent 
and there exists a peak-like component, which we model 
as a Lorentzian, near 5 x 10 -4 Hz whose strength increases 
with photon energy. 

5. CROSS-SPECTRAL PROPERTIES 

The cross-spectrum is used to quantify interband cor- 
relations as a function of temporal frequency; analogous 
to the fact that the PSD is the Fourier transform of the 
ACF, the cross-spectrum is the Fourier transform of the 
CCF. The cross-spectrum is a complex number; coherence 
and time lags can be derived from the squared magnitude 
and argument, respectively. Because these quantities in 
compact accreting objects tend to be temporal-frequency 
dependent, they cannot be studied using cross-correlation 
functions alone. Further details of the cross-spectrum and 
its relevance to XRB and AGN observations can be found 
in e.g., Vaughan & Nowak (1997), Nowak et al. (1999a), 
Papadakis, Nandra & Kazanas (2001) and VFN03. 

5.1. The Coherence Function 

The coherence function 7 2 (/) between two time series 
measures the fraction of mean-squared variability of one 
series that can be attributed to other, as a function of 
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temporal frequency. If the two time series are related by a 
simple, linear transfer function, then they will have unity 
coherence. The coherence is defined as the magnitude 
squared of the cross spectrum normalized by the product 
of each light curve’s PSD. Here, we use the discrete ver- 
sions of the cross-spectrum and coherence functions, given 
ine.g., §5.2.2 of VFN03. 

We correct the coherence function for the influence of 
Poisson noise by following Eq. 8 of Vaughan & Nowak 
(1997). We assume that the condition of high coherence 
and high powers required to use this formulation are valid. 
Specifically, the intrinsic variability power must be at least 
a factor of a few times y/m x Pp sn , where m is the coherence 
binning factor (quantified below); this condition is thus 
valid for temporal frequencies up to roughly 1 — 2 x 10 -3 
Hz. 

Coherence between the S-M, S-H and M-H band pairs 
was measured separately for each of the eight light curves, 
then combined, sorted in temporal frequency, and binned 
in groups of m=30 to form the final coherence function es- 
timate. As shown in Figure 8, the coherence is near unity 
(> 0.8 — 0.9) over the lowest temporal frequencies sampled 
for the S-M and M-H pairs. For S-H, the coherence is a 
bit lower, ~ 0.6 — 0.8, at these temporal frequencies. At 
the highest temporal frequencies sampled, the coherence 
significantly drops sharply towards zero in all cases. 

This cutoff at high temporal frequencies can be intrin- 
sic, but could also be due to Poisson noise. Monte Carlo 
simulations were used as a check on the accuracy of the 
coherence estimates and to determine if significant devia- 
tions from unity could be attributed to the effect of Poisson 
noise. For the 2005 long-look, two light curves of duration 
several Msec were simulated using the same random num- 
ber seed and assuming identical PSD shapes (the average 
of the best-fitting singly-broken PSD model for each band 
pair) to produce two light curves with intrinsic unity co- 
herence. The simulated light curves were rebinned and 
resampled to match the observed light curves (e.g, broken 
into seven smaller light curves matching the sampling of 
the observed ones), and rescaled to match the mean ob- 
served count rates in the softer and harder bands. Poisson 
noise was then added to the data according to the Pois- 
son distribution, i.e., data were randomly deviated using 
a Gaussian convolved with the square root of the observed 
mean error squared divided by the sampling time. Simu- 
lated light curves for the 2001 data were produced in the 
same way. One hundred sets of simulated data were pro- 
duced. For each set, the coherence and its uncertainty was 
calculated as for the real data (e.g., combining the eight 
simulated light curves). At most temporal frequencies be- 
low ~ 10~ 3 Hz, the simulated coherence was reasonably 
close to unity, and the uncertainty in the coherence was in 
reasonable agreement with the scatter, indicating that the 
coherence estimation was reasonably accurate. The 90% 
and 95% lower limits to the coherence were determined. 
Any measurement of the observed coherence function ly- 
ing below these lines thus indicates a drop in coherence 
that is intrinsic to the source, and not an artifact of the 
Poisson noise, at 90% or 95% confidence. For most tem- 
poral frequencies below ~ 10~ 3, ° Hz, the artificial drop in 
coherence is small. 

The energy-dependence of the high-frequency PSD slope 


could also contribute to a reduction in measured coher- 
ence. The Monte Carlo simulations were repeated, assum- 
ing the best-fit sharply- broken power-law model shapes. 
The effects of Poisson noise were again added. The sim- 
ulated coherence was typically only ~ 1%, 3 — 4% and 
1 — 3% lower for the S-M, S-H and M-H bands, respec- 
tively, compared to simulations assuming identical PSD 
shapes. The effect of having differing PSD shapes is thus 
small. The 90% and 95% lower limits to the coherence were 
determined and are plotted in Figure 8. The observed de- 
viations from unity are generally larger than the artificial 
drop due to Poisson noise. The decrease in intrinsic coher- 
ence as the energy separation of the bands increases and 
as temporal frequency increases towards ~ 10“ 3 Hz is also 
likely intrinsic. 

Having quantified the decrease in coherence e.g., due 
to Poisson Noise at each temporal frequency, we added 
the mean value of the drop in coherence from the Monte 
Carlo simulations to the observed coherence to obtain 
an estimate of the ’intrinsic’ coherence, plotted in the 
bottom panel of Figure 8. We fit the ’intrinsic’ coher- 
ence functions with exponential decay functions of the 
form 7 2 (/) = Cexp(— (/// c ) 2 ), where f c is the tempo- 
ral frequency cutoff in Hz and C is a constant, ignoring 
points above 1.2 x 10~ 3 Hz, and using least-squares fit- 
ting, assuming equal weighting to all points. We find that 
l°g(/ c ) for the S-M, S-H and M-H coherence functions are 
—2.82 ±0.03, —2.80 ±0.03 and —2.81 ±0.03, respectively. 

5.2. Temporal Frequency- Dependent Time Lags 

Cross-correlation functions (CCFs) for XRBs and AGN 
generally show a peak near zero lag but are often asym- 
metrically skewed towards hard lags (e.g., M c Hardy et 
al. 2004), suggesting there is a soft-to-hard time lag of 
some sort, but it is not a simple delay. We measured 
the CCFs for the 2005 data using e.g., the Interpolated 
Correlation Function (White & Peterson 1994). Because 
large gaps such as those between revolutions complicate 
analysis, we calculated CCFs for each individual revolu- 
tion separately. As expected, all peak lags are consistent 
with zero and CCFs exhibit a slight soft-to-hard asym- 
metric skew. The logical direction is to explore temporal 
frequency-dependent time lags via the cross-spectrum. 

Obtaining sensible measurements of time lags via the 
cross-spectrum is difficult for observations of AGN, how- 
ever; it requires the combination of a large quantity of 
uninterrupted data with a high sampling time resolution 
and a high variability-to-Poisson noise ratio. VF03 did not 
present lags for the 2001 data; we also attempted to derive 
time lags from the 2001 pn data and did not detect any 
obvious time lag trends, likely due to the combination of 
intrinsically small lag values (see below), high scatter, and 
low temporal frequency resolution. 

We used Eqs. 28 and 29 in §5.3.1 of VFN03 to calcu- 
late the phase lags 4>{f)\ time lags were then obtained as 
r(/ ) = <j>(f)/2nf. Time lags were calculated for each of 
the eight light curves separately, then combined, sorted 
in frequency, and binned in groups of 50. Errors on 
phase lags were calculated using Eq. 16 of Nowak et al. 
(1999a); in this formulation we used the coherence func- 
tion appropriate for each binned point. Sensible lags (e.g., 
lags all have the same sign and show a trend) were ob- 
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tamed for the 4 lowest temporal frequency bins, spanning 
4.4 x 10~ 5 — 3.3 x 10~ 4 Hz. The results are shown in Fig- 
ure 9. Errors are large, primarily because coherence in this 
temporal frequency range is not always exactly unity. Us- 
ing Eq. 16 in Nowak et al. (1999a), we estimated the effec- 
tive sensitivity limit on time lag detections due to Poisson 
noise, using the PSD shape and normalization from best- 
fit M-band singly-broken PSD model shape in §4. Limits 
are shown as dotted lines in Figure 9. 

Monte Carlo simulations, based on the simulated light 
curves used for the coherence functions, were performed 
to verify that the time lags were not an artifact of Poisson 
noise nor of the differing PSD slopes between two bands. 
For all simulations, time lags were <24 s for the two lowest 
temporal frequency bins and <15 s for the next two bins. 
We can reasonably conclude that the observed time lags 
are not an artifact of Poisson noise or PSD shape differ- 
ences, and are henceforth assumed intrinsic to the source. 

The time lag spectrum displays the general qualities 
seen in many previous time lag spectra: a general increase 
towards lower temporal frequencies and as the energy sep- 
aration of the bands increases. Assuming that the time 
lags are proportional to r(/) = C/ _1 for simplicity allows 
us to calculate fractional time lags, i.e., the lags’ fraction 
of the temporal frequency. We find that C = 0.002, 0.009, 
and 0.007 for S-M, S-H and M-H, respectively. These val- 
ues of C are roughly correlated with values of \og(E 2 / E\) 
= 0.32, 1.03 and 0.72, respectively, where E± and ffy are 
the relatively softer and harder average photon energies 
in each band (§4.2.1). Similar correlations have been re- 
ported previously in XRBs and in the NL Syl Ark 564 
(Arevalo et al. 2006a). The physical implications of these 
time lags and comparisons to other objects will be dis- 
cussed further in §6. 


6. DISCUSSION 

An XMM-Newton observation of Mkn 766 spanning six 
revolutions in 2005 caught the source in a relatively low 
flux state; with the average 0.2-12 keV flux about half 
of the flux during the 107 ksec XMM-Newton observation 
in 2001. Reduction of the 2005 observation yielded seven 
long, virtually uninterrupted light curves of this highly- 
variable source totaling 442 ksec of good exposure time. 

We find no evidence for the underlying variability pro- 
cess to have changed from 2001 to 2005, despite a factor of 
2 difference in average flux. Specifically, the PSD shapes 
and coherence functions between 2001 and 2005 are consis- 
tent and both observations appear to have consistent linear 
rms-flux relations and constant <F vav > values. The Mkn 
766 light curves are thus stationary over a four year dura- 
tion as far as the second-order moments are concerned. 

Combining the 2005 and 2001 observations, we have 
constructed a binned PSD with the highest temporal fre- 
quency resolution in log space for an AGN PSD to date. 
We find energy-dependent behavior not reported previ- 
ously in AGN PSDs, as discussed in §6.1. In addition, 
time lags have been detected in Mkn 766 for the first time; 
§6.2 discusses the cross-spectral results and compares them 
to those of other compact accreting objects. Finally, the 
results are discussed in the context of phenomenological 
variability models in §6.3. 


6.1. Summary of PSD Results 

The high temporal frequency resolution of the PSDs 
allows us to place tight constraints on localized features 
(a PSD break or the presence of a peak-like feature, i.e., 
a QPO) in any model we fit. Our initial approaches to 
model the PSD centered on broken power models. Model 
fits cannot formally distinguish between a sharp break or 
a slow rollover, despite the high frequency resolution; both 
models fit the data well. In general, there are no ob- 
vious trends in the power-law slope above or below the 
break, suggesting the existence of a “universal” PSD shape 
at all bands whose break scales with photon energy. In 
that case, fit parameters may depend on the temporal fre- 
quency range being studied. Our main result is that fy 
increases with photon energy Neither the increase in fy 
nor an energy-independent PSD shape has been conclu- 
sively demonstrated previously in AGN. 

6.1.1. The increase in fy with photon energy 

Our estimate for the T-band break frequency is 
log(/ b ) = —3.36 ± 0.04 (sharply-broken model, all param- 
eters free). This value is consistent with that of VE03, but 
the uncertainty is much lower compared to VF03 or even 
to VFN03 for the three XMM-Newton revolutions of data 
for MCG-6-30-15. For the first time in an AGN PSD, the 
break frequency is seen to increase with energy, a result 
that is robust to the specific broken power-law model shape 
used. The break frequency tends to increase by ~ 0.2 — 0.4 
in the log between the S and H bands (a factor of ~10 in 
average photon energy). 

Most previous Seyfert PSDs simply did not have suffi- 
cient temporal frequency resolution to detect significant 
changes in /b with photon energy such as those seen here. 
It is plausible that for a limited dynamic range PSD with 
low frequency resolution, the effect of having fy intrin- 
sically increasing with photon energy could lead instead 
to apparent measurements of the power-law slope above 
the break, f3, flattening with energy. Indeed, when we fit 
the Mkn 766 PSD with 7 and fy frozen, /3 flattens with 
energy, suggesting that previously observed trends of (3 
with photon energy may instead be due to an intrinsic en- 
ergy dependence of /b- Alternatively, it could be the case 
that this phenomenon might exist only in Mkn 766; high- 
quality, high-resolution PSDs are needed for a sample of 
Seyferts to determine whether this is the case. 

The fact that the observed increase in fy with energy is 
small suggests that PSD measurement spanning a range 
of photon energies wider than those used here would be 
needed to further separate fy in energy. In the case of 
XRBs, it is possible that this phenomenon may not have 
been reported due to insufficient spread in the hard X-ray 
photon energies studied. 

6.1.2. The possible presence of a Lorentzian component 

The alternate to modeling the PSD with /b as a free 
parameter is to assume that fy is energy-independent. 
Achieving a good fit requires the presence of a peak-like 
component, which we model as a Lorentzian with Q — 5, 
that is centered at 4.6 ± 0.4 x 10~ 4 Hz and whose rms 
variability R increases from < 1.5% in the soft band to 
5.5 ± 0.5% in the hard band. Such a detection, though 
model-dependent, is also a new result for AGN. 
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Following the definition of van der Klis (2006), the 
Lorentzian component can be identified as a QPO. Most 
previous reports of deterministic behavior or claims of 
QPOs in AGN have been refuted, have failed the test of 
repeated observations, or did not use the appropriate sta- 
tistical significance tests (see Vaughan & Uttley 2005 and 
Vaughan 2005). We our basing our claim on modeling of a 
PSD that is properly binned, has well-defined errors, has 
very high temporal frequency resolution, and is derived 
from a large quantity of high-quality data. We sample 
almost 300 cycles of this feature throughout the duration. 

We reiterate that our claim of the presence of a QPO in 
the PSD is highly model-dependent, since a simple broken 
power-law model can fit the data well. The requirement 
for the QPO arises only when /b is assumed to be energy- 
independent. Both models adequately describe the energy 
dependence of the PSD well, though Occam’s Razor may 
have one prefer the simple, broken-power law model, since 
it consists of only the broadband PSD component and does 
not make any a priori assumption about the behavior of 
/b with energy. 

Assuming that the QPO at 4.6 x 10~ 4 Hz is real, there 
are intriguing similarities to those observed in XRBs. 
QPOs whose rms variability R increases with photon en- 
ergy have been observed in XRBs. For example, the 450 
Hz QPO in GRO J1655-40 and the 40 Hz QPO in GRS 
1915+105 are detected only above 13 keV by Strohmayer 
et al. (2001a) and (2001b), respectively. In Mkn 766, the 
value of R is highly model-dependent, but also increases 
with photon energy and displays roughly similar values, 
from < 1.5% at ~0.4 keV to 5.5% at ~4 keV. It is not un- 
reasonable to compare the (relatively softer) energy bands 
in Mkn 766 with the (harder) energy bands in XRBs, since 
softer bands in XRBs may be dominated by thermal disk 
emission. 

Photon indices in the low/hard, high/soft and very high 
state in XRBs tend to be ~ 1.5 — 2.0, >2 — 3, and 
~ 2.4 — 2.9, respectively (e.g., McClintock & Remillard 
2006; Homan et al. 2005). The photon index of the power- 
law continuum in Mkn 766 is ~ 2.3 — 2.4 (e.g., Turner 
et al. 2006, Miller et al. 2006). A direct comparison be- 
tween Mkn 766 and XRBs in the very high state therefore 
might not be rejected on those grounds. QPOs, however, 
tend to appear only in XRB states when the thermal disk 
emission does not dominate the energy spectrum (van der 
Klis 2006, McClintock & Remillard 2006). The existence 
of a QPO in Mkn 766 would therefore argue against Mkn 
766 being an analog of the (thermal-dominated) high/soft 
state in XRBs. 

One may ask whether the possible QPO in Mkn 766 is 
the analog of low-frequency (LF) or high-frequency (HF) 
QPOs. Type A and B LF QPOs usually appear at 4-9 Hz. 
Assuming that QPO frequencies scale linearly with black 
hole mass, and assuming 10 M© for the mass of stellar- 
mass black holes, the mass of the black hole in Mkn 766 
is thus estimated to be 1 — 2 x 10 5 M©. This is somewhat 
low compared to estimates of Mkn 766’s black hole mass: 
Botte et al. (2005) measure a stellar velocity dispersion <r* 
of 81 ± 17 km s -1 in Mkn 766. Using the relation between 
u* and Mbh as parameterized by Tremaine et al. (2002), 
this suggests Mbh = 3. 6^9 x 10 6 M©. Estimates using 
the width of the H/J line and the Kaspi et al. (2000) re- 


lation lie in the range 0.8 — 10.0 x 10 6 M© (Wandel 2002; 
Bian & Zhao, 2003a, 2003b, 2004; Botte et al. 2005). On 
the other hand, that fact that XRB QPOs tend to appear 
only in states when the thermal disk emission does not 
dominate the energy spectrum suggests that QPOs are as- 
sociated with the corona instead of the disk (Reig et al. 
2006; van der Klis 2006, McClintock & Remillard 2006). 
QPO frequencies may thus scale with some other prop- 
erty besides or in addition to black hole mass. Type C LF 
QPOs tend to have R values which may be too high to 
compare to Mkn 766. 

HF QPOs, usually are seen near 100-300 Hz, and usu- 
ally with R ~ 1 — 3%. They are often found in pairs, at 
frequencies of 2 and 3 times an unseen “fundamental” fre- 
quency /o- The Mkn 766 PSD does not have the resolution 
to determine if the QPO corresponds to the 2/o or the 3/o 
peak being dominant, or if it is a blend of the two. Based 
on observations of three black hole XRBs, McClintock & 
Remillard (2006) have suggested that /o scales linearly 
with black hole mass Mbh as 931 Hz (Mbh/M©) _1 , as- 
suming similar values of the spin parameter; this linear 
scaling suggests that HF QPOs in black holes may repre- 
sent disk oscillations described by general relativity (see 
Abramowicz & Kluzniak 2001). Extrapolating this rela- 
tion to Mkn 766 and assuming /o = 1.8 — 2.4 x 10~ 4 Hz, 
this implies a black hole mass of 3.9 — 5.2 x 10 6 M©, con- 
sistent with the optically-estimated black hole masses. On 
the other hand, HF QPOs are not always found very close 
to the break in the broadband component, as is the case 
here. 

In the end, the community needs very high-quality data 
on a reasonably sized sampled of AGN to be able to con- 
firm or refute the presence of QPOs in AGN in order to 
shore up the comparisons between AGN and XRBs. The 
Mkn 766 observations are an important step towards this 
goal, demonstrating that high-quality data in the form of 
a collection of long, uninterrupted exposures can success- 
fully yield an enormous payoff for AGN variability inves- 
tigations. 

6.2. Summary of cross-spectral results and comparison to 

other objects 

The coherence is generally flat at the lowest temporal 
frequencies probed, with values typically ~ 0.8 — 0.9 for 
S-M and M-H and ~ 0.6 — 0.8 for S-H. The coherence 
drops towards zero near ~ 10~ 3 Hz; Monte Carlo simula- 
tions prove that these deviations from unity are intrinsic 
to the source and not the result of Poisson noise or energy- 
dependent PSD shapes. The decrease in coherence as en- 
ergy separation increases has been reported previously for 
Mkn 766 by VF03; we in fact see values of the coherence 
similar to VF03. This similarity again supports the con- 
clusion of §3, that there is no evidence for the underlying 
variability process to have changed significantly between 
2001 and 2005 despite the overall lower average flux in 
2005. 

In other black hole systems, the coherence at tempo- 
ral frequencies below the break is generally between 0.8 
and 1.0, though coherence between bands similar to the S 
and H bands used here can be as low as ~ 0.6 — 0.7 (e.g., 
VFN03 & M c Hardy et al. 2004); Mkn 766’s coherence func- 
tions thus agree with those in other compact objects. As 
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explained in Vaughan & Nowak (1997), coherence between 
two energy bands can be lost if the transfer function relat- 
ing the light curves is non-linear, or if there are multiple 
uncorrelated flaring regions. The fact that coherence in 
Mkn 766 and several other compact objects is not exactly 
unity may support some degree of non-linearity present in 
the transfer functions relating relatively softer and harder 
bands, and which increases for larger energy separations. 

The coherence function between the S and H bands 
drops to ~l/e around 1 x 10~ 3 Hz in both MCG-6-30- 
15 (VFN03) and Mkn 766, while the drop in NGC 4051’s 
S-H coherence occurs at roughly ~ 4 x 10~ 3 Hz (M c Hardy 
et al. 2004). MCG-6-30-15’s black hole mass is roughly 
3 x 10 6 M© (M c Hardy et al. 2005), similar to that for Mkn 
766, while that for NGC 4051 is 1.9 ± 0.8 x 10 6 M© (Pe- 
terson et al. 2004). A larger sample of coherence func- 
tions of Seyferts spanning a larger black hole mass range 
is needed to verify if coherence time scales scale with black 
hole mass, just as PSD break time scales do. 

As explained by e.g., Nowak & Vaughan (1996) and 
Nowak et al. (1999a), the sharp drop in coherence at tem- 
poral frequencies above the PSD break could indicate that 
the corona is dynamic on such short time scales. Such 
short time scales could be associated with formation time 
scales or changes in temperature and/or the physical struc- 
ture of the corona. The relatively high levels of coherence 
at low temporal frequencies could indicate that the corona 
is effectively static on long time scales. 

We measure temporal-frequency dependent phase lags 
in Mkn 766 for the first time; lags increase towards lower 
temp frequency and as energy separation of bands in- 
creases, as seen in other black hole systems. Fractional 
time lags are about 0.01-0.02 the Fourier time scale. These 
values are consistent with most values measured so far for 
other AGN and many XRBs (see e.g., Figure 8 of Arevalo 
et al. 2006a), with fractional time lags usually measured 
to be < 0.2, though Ark 564’s time lag spectrum (Arevalo 
et al. 2006a) seems to be an exception among AGN. 

Assuming a PSD break frequency of 10~ 3 ' 4 Hz, the 
Mkn 766 lag spectrum is probing temporal frequencies 
0.14 — 1.00 times the break frequency. We can compare 
Mkn 766’s fractional lags to those measured for Cyg X- 
1 within ± a decade of temporal frequency of the PSD 
breaks. In the low/hard and high/soft state of Cyg X-l, 
fractional time lags tend towards 0.03-0.1 (e.g., Nowak et 
al. 1999a), while in the intermediate state, fractional time 
lags are higher, ~ 0.1 -0.2 (e.g., Pottschmidt et al. 2000). 
On the basis of these fractional lag values, Mkn 766 is likely 
not a supermassive analog of Cyg X-l in the intermediate- 
state, though we cannot determine on the basis of lags 
alone if analogy with the low/hard or high/soft state is 
appropriate. Furthermore, the time lag spectrum is not 
of high enough quality to distinguish between an intrin- 
sic lag spectrum which is step-like, as in the intermediate 
state of Cyg X-l, or adheres to a form more resembling a 
power-law. 

6.3. Phenomenological Variability Models 

Many models have been invoked to explain the red 
noise variability properties of Seyferts and XRBs, includ- 
ing shot-noise models (e.g., Merloni & Fabian 2001), ro- 
tating hot-spots on the surface of the accretion disk (Bao 


& Abramowicz 1996), self-organized criticality (“pulse 
avalanche” models; Mineshige, Ouchi & Nishimori 1994), 
magnetohydrodynamical instabilities in the disk (Hawley 
& Krolik 2001), and inwardly- propagating fluctuations in 
the local accretion rate (Lyubarskii 1997). As noted by Ut- 
tley, M c Hardy & Vaughan (2005), shot-noise models have 
difficulty reproducing the linear rms-flux relation, and, as 
noted by Vaughan & Nowak (1997), they lead to low ob- 
served coherence values unless each separate emitting re- 
gion independently has the same linear transfer function 
between energy bands. 

Thermal Comptonization of soft seed photons by a hot 
corona is a likely explanation for the X-ray emission. It is 
supported by X-ray energy spectra, and the soft-to-hard 
time lags could be attributed to the time scale for seed 
photons to diffuse through the corona and undergo multi- 
ple up-scatterings. However, the simplest Comptonization 
models do not predict time lags that depend on temporal 
frequency, as observed. 

A model in which inwardly-propagating variations in 
the local mass accretion rate rh are responsible for the 
observed X-ray variability seems to be able to explain 
many of the observational results in Seyferts and XRBs 
(Lyubarskii 1997, Churazov, Gilfanov & Revnivtsev 2001), 
Kotov et al. 2001, Arevalo & Uttley 2006; see also VFN03 
and M c Hardy et al. 2004). In this model, variations at 
a given radius are associated with the local viscous time 
scale, so relatively smaller radii are associated with rela- 
tively more rapid variations. Perturbations in m propa- 
gate inward, until they reach, and modify the emission of, 
the central X-ray emitting region. The resulting PSD has 
a 1// form across a broad range of temporal frequencies, 
the net sum of variations from a wide range of radii. How- 
ever, if the X-ray emitting region is extended, it essentially 
acts like a low-pass filter on these 1// variations: with rel- 
atively higher-frequency variations imprinted on the light 
curve originating only from the smallest radii, this pro- 
duces a cutoff in the high-frequency PSD, as is commonly 
observed. Kotov et al. (2001) suggested that the spec- 
trum can be a function of radius, with relatively harder 
emission emanating from smaller radii, and therefore as- 
sociated with more rapid variability. This class of models 
thus predicts several additional properties which are ob- 
served in Seyferts and XRBs, namely the presence of the 
linear rms-flux relation, and the dependence of coherence 
and time lags on temporal frequency and on the energy 
separation of the bands. 

Arevalo & Uttley (2006) show that if relatively harder 
X-ray bands are associated with emissivity profiles that 
are more centrally concentrated, and if the PSD break is 
at least partly due to filtering produced by having an ex- 
tended X-ray emission region, then the harder bands will 
retain more variability power at the highest temporal fre- 
quencies, leading to break frequencies in the broadband 
PSD which increase with increasing photon energy. We 
note that the PSDs of Mkn 766 as shown in Figure 5b, 
bear a striking qualitative resemblance to the simulated 
PSDs in Figure 2 of Arevalo & Uttley (2006). The de- 
gree of curvature in the Mkn 766 PSDs is uncertain, and 
rollover frequencies are relatively close together, so it is 
not straightforward to relate break frequencies to emissiv- 
ity index or profile. However, we note that the ’soft’ and 
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’hard’ PSDs in Figure 2 of Arevalo & Uttley (2006), sim- 
ulated assuming emissivity indices which differ by 2, yield 
break frequencies which are separated by about a decade. 
In Mkn 766, the S- and H-band break frequencies are not 
as separated, suggesting that, in the context of this model, 
the emissivity indices for the S- and H- bands differ by less 
than 2. The relatively small fractional time lags observed 
in Mkn 766 (0.01-0.02 for the S-H bands) are consistent 
with a difference in emissivity indices of ~0.5 to a few. 

7. CONCLUSIONS 

We have analyzed the energy-dependent variability 
properties of the Narrow Line Seyfert 1 Mkn 766 using 
light curves obtained from two XMM-Newton observa- 
tions. An observation spanning six revolutions in 2005 
yielded seven uninterrupted EPIC light curves totaling 442 
ksec of good exposure time; an observation in 2001 yielded 
an eighth uninterrupted light curve with an exposure of 
107 ksec. A linear rms-flux relation, or equivalently, con- 
stant <i ? var >, was reported for the 2001 observation by 
Vaughan et al. (2003a). These relations are confirmed for 
the 2005 observation, and shown to be consistent with the 
corresponding relations from the 2001 observation. Com- 
paring the 2001 and 2005 PSDs shows the PSDs to be con- 
sistent within the expected scatter. Variability properties 
between the 2005 and 2001 observations are thus consis- 
tent with each other, despite a 0.2-12 keV flux in 2005 
which is about half that of the 2001 observation. We com- 
bined the 2001 and 2005 data sets in order to probe power- 
and cross-spectral properties with temporal frequency res- 
olution that is unprecedented for an AGN. 

We have explored two parameterizations of the energy- 
dependent PSD. Both fit the data equally well, but in ei- 
ther case, they signify behavior never seen previously in an 
AGN PSD. Fitting simple broken power-laws, we find PSD 
breaks at nearly the same temporal frequencies as VF03, 
but the break frequencies /b are seen to increase with pho- 
ton energy, typically by ~0.2-0.3 in the log between ~0.4 
keV and ^4 keV emission. There is no obvious trend in 
PSD power-law slope, and the data are consistent with a 
universal PSD shape whose break frequency increases with 
photon energy. Previous Seyfert PSDs lacked the tempo- 


ral frequency resolution to cleanly detect the energy de- 
pendence in /b- Keeping /b fixed with energy does not 
provide a good description of the energy-dependent PSD. 
However, we can achieve a good fit using a sharply-broken 
power-law, where the break frequency and power-law slope 
above the break are roughly energy-independent, but we 
must include a narrow peak component, i.e., a QPO, at 
4.6 ± 0.4 x 10~ 4 Hz. Interestingly, the rms variability am- 
plitude of this QPO increases with photon energy in a 
manner similar to QPOs seen in the PSDs of black hole 
XRBs. 

Cross-spectral properties are qualitatively consistent 
with those measured for black hole systems previously. 
The coherence is generally flat at the lowest temporal fre- 
quencies probed, with values typically ~ 0.6 — 0.9. The co- 
herence drops towards zero at temporal frequencies higher 
than roughly 10~" 2 ' 8 Hz; Monte Carlo simulations indicate 
that the drop is intrinsic to the source and not an artifact 
of Poisson noise or differing PSD shapes. 

Temporal frequency-dependent time lags are measured 
from the cross-spectrum for the first time in Mkn 766. As 
is the case with time lags in other objects, lags increase 
as temporal frequency increase and as energy separation 
of the light curves increases. Fractional time lags are typ- 
ically <0.01, consistent with many other Seyferts as well 
as with the high/soft and low/hard states of Cyg X-l. 

The results were discussed in the context of several vari- 
ability models, including models incorporating inwardly- 
propagating fluctuations in the local mass accretion rate, 
e.g., Lyubarskii (1997). Notably, the observed increase 
in the PSD break frequency with photon energy is qual- 
itatively consistent with the prediction put forth by e.g., 
Arevalo & Uttley (2006), who model variations in a disk 
with relatively harder X-ray emission being more centrally 
concentrated. 


We thank the referee for providing useful comments. 
This work has made use of observations obtained with 
XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA member states 
and the US (NASA). 


APPENDIX 
PSD CONSISTENCY 

We tested for strong non-stationarity in the form of a change in the shape of the underlying PSD, between the 2001 and 
2005 PSDs and between revolutions in the 2005 long-look. For two PSDs observing the same stationary process at two 
different times, the only difference in PSD values at a given temporal frequency should be consistent with the expected 
scatter. The scatter in the measured PSD about the intrinsic, underlying PSD follows a y 2 distribution with two degrees of 
freedom (e.g., Priestley 1981); Papadakis & Lawrence (1995) thus outlined a method to test for non-stationarity, defining 
a statistic S to quantify differences between two PSDs. For a stationary process, S will be distributed with a mean of 0 
and a variance equal to 1 (see also VFN03). 

We considered PSDs for each revolution from 999 to 1003A using the summed pn + MOS1+MOS2 T, S, M and H-band 
light curves as well as for the summed pn + MOS2 light curve from 2001. Light curves were binned to 60 seconds and 
truncated to a common duration of 52.6 ksec; PSDs were generated by binning periodogram points by a factor of 25. 

The S statistic was calculated by summing over temporal frequencies <10 -2 ' 6 Hz to exclude PSD points that were 
dominated by power due to Poisson noise. For all PSD pairs and bands, |Sj was found to 0.54 or less. Extending the 
duration to 74.0 ksec to compare revolutions 265 and 999-1002 and extending the duration to 92.2 ksec to compare 
revolutions 265 and 1000-2 also yielded similarfy small values of (FI (<0.5) for all bands. To summarize, all values of |Sj 
were less than 1 for all light curve pairs and bands; we therefore have found no evidence at la confidence or greater for 
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strong non-stationarity in the PSDs of Mkn 766 at any band from one revolution to the next in the 2005 long-look, or 
between 2001 and 2005. 
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Table 1 

Unbroken Power-law Model Fits to the PSD 


Band 

X 2 /dof 

0 

A 

T 

299.3/180 

2.34±0.06 

0,99±0.03 

S 

137.4/180 

2.14±0.08 

0.79T0.04 

M 

258.3/180 

2.22±0.06 

1.00±0.03 

H 

278.4/180 

2.14±0.08 

1.16±0.03 


Note. — A is the PSD amplitude in units of Hz 1 at fo = 10 350 


Hz. 
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Table 2 


Sharply- and Slowly Broken Power-law Model Fits to the PSDs 


Band 

X 2 / do f 

7 

l°g(/b(Hz)) 

P 

A 

Sharply-Broken Power-Law Model Fits 

T 

207.5/178 

1.48i0.14 

-3.36i0.04 

3.00±0.12 

0.95±0.03 

S 

217.9/178 

1.56i0.18 

-3.55i0.09 


1.08±0.04 

M 

173.0/178 

1.54i0.14 

-3.34i0.05 

3.10±0.18 

0.91i0.03 

H 

178.0/178 

1.28i0.12 

-3.29i0.04 

3.02i0.16 

1. 01±0. 03 

T 

208.6/180 

1.46 (fixed) 

-3.39±0.05 

2.89 (fixed) 

0.99i0.03 

S 

222.6/180 

1.46 (fixed) 

-3.48i0.06 

2.89 (fixed) 

1.00i0.04 

M 

175.1/180 

1.46 (fixed) 

-3.40i0.05 

2.89 (fixed) 

1.03i0.03 

H 

187.2/180 

1.46 (fixed) 

-3.27±0.05 

2.89 (fixed) 

0.92i0.03 

T 

207.6/180 

1.46 (fixed) 

-3.36 (fixed) 

3.00i0.14 

0.94i0.03 

S 

231.6/180 

1.46 (fixed) 

-3.36 (fixed) 

2.94i0.24 

0.73i0.04 

M 

173.4/180 

1.46 (fixed) 

-3.36 (fixed) 

3.08l°; 22 

0.97i0.03 

H 

192.5/180 

1.46 (fixed) 

-3.36 (fixed) 

2.76±0.14 

1.10i0.03 

Slowly-Broken Power-Law Model Fits 

T 

210.8/178 

1.26i0.12 

-3.34±0.05 

3.34i0.12 

-3.10i0.03 

S 

218.3/178 

1.04±0.24 

-3.77i0.09 

2 . 68 ± 0.12 

-2.16i0.04 

M 

172.9/178 

1.42±0.12 

-3.23i0.06 

3.72±0.22 

-3.73i0.03 

H 

186.4/178 

1.16±0.10 

-3.20±0.04 

3.58i0.20 

-2.70i0.03 

T 

210.9/180 

1.21 (fixed) 

-3.36±0.05 

3.33 (fixed) 

-2.91i0.04 

S 

225.5/180 

1.21 (fixed) 

-3.44i0.07 

3.33 (fixed) 

-3.03i0.03 

M 

174.5/180 

1.21 (fixed) 

-3.39±0.05 

3.33 (fixed) 

-2.86i0.03 

H 

188.3/180 

1.21 (fixed) 

-3.23i0.04 

3.33 (fixed) 

-2.85i0.03 

T 

211.0/180 

1.21 (fixed) 

-3.19 (fixed) 

3.36i0.14 

-2.92i0.03 

S 

229.6/180 

1.21 (fixed) 

-3.19 (fixed) 

3.30i0.26 

-3.13i0.04 

M 

174.3/180 

1.21 (fixed) 

-3.19 (fixed) 

3.46±0.20 

-2.90i0.03 

H 

195.3/180 

1.21 (fixed) 

-3.19 (fixed) 

3.08±0.16 

-2.76i0.03 

Note. - 

— A is the PSD amplitude 

in units of Hz -1 

at the break frequency fb- 



Table 3 

PSD Model Fits Using a Sharply-Broken Power-law plus a 

Lorentzian 

Band 

X 2 /dof 

7 

log(/ b (Hz)) 

P 

A 

l°g(/L) 

R (%) 

T 

S 

M 

H 

212.7/178 

218.2/178 

178.1/178 

193.7/178 

1.46 (fixed) 
1.46 (fixed) 
1.46 (fixed) 
1.46 (fixed) 

-3.60l° O6 

-3.60l°. os 

-3-60l° os 

-3-60l°. os 

2 . 68 i 0.10 

2.5410.14 

2.7010.14 
2.42i0.12 

1.3310.03 

1.1910.04 

1.3510.03 

1.3810.03 

-3.34i0.04 

-3.34i0.04 

-3.34i0.04 

-3.34i0.04 

3.010.5 

< 2.0 
0 r+O.S 
J,o - 1.0 

5.010.5 


Note. — Results to fitting the PSDs with a model consisting of a sharply-broken power-law plus a Lorentzian. A is the PSD 
amplitude in units of Hz -1 at the break frequency fb ■ The break frequency fb was constrained to be 10 ~ 3 - 60 Hz or lower. The 
low-frequency slope 7 was fixed at -1.46. Based on the sum of the x 2 values for the S, M and H bands, the best-fit value of the 
Lorentzian centroid frequency ft was found to be io ~ 3 - 34±0 - 04 Hz; fit parameters above correspond to that value of ft- In all 
fits, the quality factor of the Lorentzian Q was assumed to be 5. 
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Fig. Al. — 0.2-12 keV EPIC light curves, summed over the pn and MOS cameras, binned to 60 s, for the 2005 observation. All panels have 
the same scale. 



Observed Energy (keV) Observed Energy (keV) Observed Energy (keV) 

Fig. A2. — F V ar spectra, constructed using EPIC-pn data, for the entire 2005 observation and for the lowest- and highest-flux revolutions. 
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Flux (0.2-1 2 keV ct S" 1 ) 

Fig. A3. — The 2005 pn data (black circles) and the 2001 pn data (gray open squares) are both consistent with a linear relation between 
rms variability amplitude (top) and with constant average fractional variability amplitude <F var >. The black solid and gray dashed lines are 
the best-fit relations for the 2005 and 2001 data, respectively. 
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10- 4 10~ 3 

Temporal Frequency (Hz) 


Fig. A4. — a) Power density spectrum for the T-band light curve. The solid line denotes the best-fit sharply-broken power-law fit. b) 
Residuals to the best-fit unbroken power-law model; residuals are measured in log space, e.g., based on the logs of the PSD model, data and 
errors, c) Residuals to the best-fit sharply- broken power-law model, d) Residuals to the best-fit slowly-bending power-law model. Power due 
to Poisson noise was included in the models in all panels. 
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Temporal frequency (Hz) 

FlG. A5. — a) Power density spectra for the S, M, and H-band light curves, plotted in / — / x P(f) space. The solid lines show the best-fit 
sharply-broken power-law model assuming the power-law slopes above and below the break /? and 7 are equal in all bands. The models 
are plotted with the constant level of power due to Poisson noise omitted for clarity, b) Same as a), but the solid lines show the best-fit 
slowly-bending power-law model assuming the power-law slopes above and below the break 0 and 7 are equal in all bands, c) Residuals to the 
best-fit unbroken power-law model, d) Residuals to the best-fit sharply-broken power-law model, e) Residuals to the best-fit slowly-bending 
power-law model (power due to Poisson noise was included in the models in panels b)-d). 
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10- 4 10- 3 

Temporal Frequency (Hz) 


Fig. A6. — a) Power density spectra for the S, M, and H-band light curves, plotted in / — fxP(f) space. The solid lines show the best-fit 
model consisting of a sharply-broken power-law, with break frequency constrained to be Hz or lower, plus a Lorentzian with a centroid 

frequency of 10 -3 - 32 Hz and quality factor Q = 5. The models are plotted with the constant level of power due to Poisson noise omitted for 
clarity, b) Residuals to the best-fit model (with power due to Poisson noise included). 



- 3.45 - 3.4 - 3.35 - 3.3 - 3.25 


i°g(f L (Hz)) 

Fig. A7. — Derived confidence contours of the Lorentzian amplitude R versus centroid frequency /l for each band; black, red, green, 
and blue denote T, S, M and H bands, respectively. 68, 90, 95 and 99% contour levels are shown (solid, dashed, dash-dotted, and dotted, 
respectively). The blue cross marks the best-fit value for the H band. The green cross marks the best-fit value for both the M and T bands. 
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5x1 0 -5 10~ 4 2x10“* 5x1 0~ 4 10~ 3 5x1 0’ 5 1 0' 4 2x10-“ 5x10“ 4 10‘ 3 5x1 O' 5 1 CP 4 2x1 0- 4 5x10“ 4 10~ 3 


Temp. Freq. (Hz) Temp. Freq. (Hz) Temp. Freq. (Hz) 

Fig. A8. — The top row shows the observed coherence functions for S-M, S-H and M-H light curve pairs. The thin solid line denotes unity 
coherence. The dashed and dotted lines denote, respectively, the 90% and 95% confidence limits for spurious lack of coherence due to the 
combination of Poisson noise and the effect of differing PSD shapes as determined by Monte Carlo simulations. The bottom row shows the 
estimate of the “intrinsic” coherence (see text for details). 



Fig. A9. — Time lag spectra for Mkn 766. Positive lags denote soft band variations leading those in the hard band. A negative lag 
was detected in the second-highest temporal frequency bin in the S-M time lag spectrum, and its magnitude is denoted by an open circle. 
The dashed lines denotes the sensitivity limits for time lag detection; see text for details. The solid lines denote the best-fit relations to 

r(f) = Cf~\ 





